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SECTION  1 


INTRODUCTION 


This  report  presents  the  results  of  a  study  program  to 
investigate  concepts  for  coherent  coupling  of  multiple  GaAlAs 
diode  gain  elements  to  produce  a  laser  system  with  the  potential 
for  generating  high  average  output  power.  It  is  desired  that 
the  resulting  laser  source  have  sufficient  brightness  to  ensure 
efficient  nonlinear  wavelength  conversion,  such  as  frequency¬ 
doubling  into  the  blue-green  spectral  region. 

The  concept  of  combining  multiple  gain  modules  into  a 
single  coherent  output  beam  has  traditionally  been  approached  by 
techniques  such  as  injection  locking,  diffraction  coupling,  or  a 
combination  of  these.  By  constrast,  this  report  assumes  a  new 
approach  to  phase-locking  modular  gain  media  utilizing  the  tech¬ 
nology  area  of  nonlinear  optical  phase  conjugation.  To  under¬ 
stand  how  optical  phase  conjugation  techniques  can  be  used  for 
this  purpose,  consider  a  resonator  structure  which  employs  a 
phase  conjugate  mirror  (PCM)  as  one  element  of  the  cavity.  Such 
a  structure  is  known  as  a  phase  conjugate  resonator  (PCR). 

Figure  1  shows  a  PCR  structure  applied  to  locking  three  gain 
media.  A  generic  PCR  has  phase  and  frequency  properties  that 
are  substantially  different  from  a  conventional  resonator. 

First,  the  output  phase  profile  is  relatively  insensitive  to 
aberrations  within  the  gain  medium.1-4  In  a  conventional  reso¬ 
nator  even  small  internal  aberrations  cause  large  diffractive 
losses  and  greatly  degrade  the  output  beam.  Second,  since  we 
are  considering  PCRs  that  use  four-wave  mixing  conjugators, the 
output  frequency  is  controlled  by  the  conjugator  pumps. 

Stripped  of  details,  the  four-wave  mixer  conjugator,  through  the 
action  of  pumps  1  and  2,  couples  the  phase  and  frequency  of  the 
gain  modules  to  that  of  the  pumps.  The  phase  at  the  output 
coupling  mirror  becomes  uniform  and  results  in  a  far-field  radi¬ 
ation  pattern  possessing  a  central  lobe  and  a  specific  side-lobe 
level  as  in  any  thinned  array. 
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Figure  1 .  Generic  PCR  structure  applied  to  coupling 
modular  lasers  by  four-wave  mixing. 
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The  actual  experimental  demonstration  of  such  a  concept  is 
beyond  the  scope  of  the  present  program.  Instead,  effort  has 
been  directed  toward  several  specific  technical  areas  related  to 
coupling  multiple  diode  gain  elements,  as  well  as  some  specific 
issues  relating  to  achieving  phase  conjugation  and  efficient 
nonlinear  wavelength  conversion  within  the  800-900  nm  spectral 
range  accessible  to  GaAlAs  laser  devices.  To  enhance  the  ulti¬ 
mate  applicability  of  this  program  to  systems  applications,  it 
has  been  implicitly  assumed  that  many  diode  gain  elements  are 
fabricated  on  a  single,  practical,  monolithic  chip.  Chip  size 
or  the  total  number  of  chips  could,  in  principle,  be  scaled 
according  to  the  system  power  requirement.  Section  2  presents 
calculations  showing  the  effects  on  the  far-field  intensity 
distribution  of  either  systematic  or  random  phase  differences 
between  adjacent  gain  elements;  these  results  dramatically 
illustrate  the  beneficial  impact  of  exact  phase  locking.  Next, 
the  use  of  phase  conjugation  to  phase-lock  the  individual  emit¬ 
ters  is  presented.  Computer  calculations  showing  the  beneficial 
effects  on  the  far-field  intensity  distribution  are  reviewed. 
Next,  the  normal  modes  of  a  multiple-element  diode  chip  are 
derived.  This  information  allows  one  to  determine  phase  rela¬ 
tionships  characteristic  of  each  mode.  These  relationships  can 
be  used  to  identify  approaches  that  might  be  employed  to  favor 
one  particular  mode  in  a  practical  device.  Section  3  presents 
calculations  showing  that  a  phase  conjugate  mirror  based  on 
stimulated  Brillouin  scattering  (SBS)  in  either  a  CS2-filled 
hollow  light  pipe,  or  a  solid  fused  silica  waveguide  is  possible 
in  this  spectral  range  at  power  levels  of  only  <  1  watt.  The 
advantage  of  the  SBS  phase  conjugate  mirror  is  that  no  external 
pumps  are  required,  leading  to  an  improved  overall  system  effi¬ 
ciency.  Section  4  addresses  issues  relating  to  achieving  effi¬ 
cient  nonlinear  wavelength  conversion  with  a  diode-based  source, 
and  presents  an  approach  for  minimizing  thermal  problems  in  a 
high  average  power  frequency  doubler.  Previous  experimental 
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demonstrations  have  established  the  feasibility  of  frequency¬ 
doubling  the  output  of  diode  lasers  and  have  indicated  the 
relevant  scaling  laws.  Using  this  information,  we  estimate  that 
an  extracavity  doubling  efficiency  of  40  percent  is  possible 
using  thirteen  commerical  GaAlAs  diodes  with  non-crit ically 
phase  matched  KNb03 .  This  assumes  that  the  beam  compacting 
problem  has  been  solved. 
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SECTION  2 


COUPLED  GAIN  MEDIA 


Many  techniques  have  been  suggested  for  combining  multiple 
gain  elements  into  a  single  coherent  beam.  These  techniques 
include  injection  locking  and  path  synchronization,  diffraction 
coupling,  and  combinations  of  these.  Figure  2  illustrates  two 
systems  that  use  injection  locking  to  couple  multiple  elements. 
Figure  3  shows  a  schematic  of  diffraction  coupling  with  a 
reflective  spatial  filter.  In  general,  injection  locking  of 
multiple  diode  elements  appears  unattractive  because  of  the 
requirement  that  the  effective  path  length  of  each  element  (the 
product  of  the  refractive  index  and  the  geometric  length)  must 
be  controlled  to  within  a  small  fraction  of  a  wavelength.  Dif¬ 
fraction  coupling  has  been  investigated  by  others  5,6  and  is 
not  the  subject  of  the  current  program.  Recently,  phase  conju¬ 
gate  reflection  by  four-wave-mixing  has  developed  to  the  point 
where  it  may  be  used  to  build  a  multiple  gain-module  laser  with 
phase  and  frequency  locking.  Figure  1  showed  a  system  with  a 
four-wave-mixing  conjugator  used  to  couple  the  modules.  The 
analysis  presented  in  this  section  was  undertaken  to  provide  a 
framework  within  which  coupling  techniques  can  be  compared  and 
evaluated.  We  begin  with  a  discussion  of  the  far-field  inten¬ 
sity  distribution  that  arises  by  coupling  multiple  gain  ele¬ 
ments,  and  show  the  degradation  in  the  far— field  intensity 
distribution  that  arises  from  uncompensated  phase  differences 
between  individual  emitters.  The  next  section  illustrates  how 
phase  conjugation  compensates  such  phase  differences  to  improve 
beam  quality.  Far-field  effects  arising  from  physical  gaps 
between  emitters  (to  allow  for  heat  sinks,  for  example)  are 
presented  next,  along  with  some  basic  approaches  to  minimize 
these  deleterious  effects.  We  conclude  with  a  calculation  of 
the  normal  modes  of  a  multiple  planar  waveguide,  representing 
the  multi-element  diode  chip  assumed  to  be  the  basic  building 
block  of  these  diode  coupling  concepts. 
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b.  INJECTION  LOCKING  VIA  MUTUAL  COUPLING. 

THE  CROSS  COUPLING  ILLUSTRATED  HERE  IS  SCHEMATIC 


Figure  2.  Two  classes  of  injection  locking  techniques. 

Neither  of  these  techniques  ensures  a  co¬ 
phased  output  for  wavefronts  0^  O  ,  and  0  . 


16 


11503-7 


LASER  OUTPUT 
FIELD  DISTRIBUTIONS 


(MAGNIFIED) 

Figure  3.  Phase-locking  of  three  lasers  by 
diffraction  coupling  with  spatial 
filters. 
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2.1  FAR-FIELD  INTENSITY  DISTRIBUTION  OF  A  PHASED  LASER  ARRAY 

The  understanding  of  a  phased  laser  array  entails  the 
knowledge  of  the  relative  phases  between  adjacent  lasers.  In 
general,  it  is  a  quite  complicated  problem  due  to  the  uncer¬ 
tainty  in  the  semiconductor  laser  parameters,  i.e.,  current 
density,  index  of  refraction,  gain  homogeneity,  etc.  However, 
it  is  possible  to  estimate  the  magnitude  of  the  phase  shifts  by 
examining  the  changes  arising  from  temperature  variation,  doping 
concentration  and  gain  length.  The  optical  phase  <j>  is  defined 
as 


♦ 


(D 


where  XQ  is  the  vacuum  wavelength,  n  is  the  index  of  refraction 
and  L  is  the  propagation  distance.  The  change  in  <f>  is  then 


or. 


A  4> 


2ir  r  _ 
=  ‘ —  (A  nL 

A  o 


+  nAL) 


(2) 


A<t>  =  |[  (A  n)  t  +  (An)c]  L  +  nAL} 


(3) 


where  (in)t  and  (An)c  are  the  contributions  due  to  the  tempera¬ 
ture  and  dopant  concentration  respectively.  In  a  semiconductor 
laser,  the  contribution  to  An  due  to  a  temperature  change  is 
estimated7  to  be  dn/dT  =  3  x  10~V°K  at  850  nm.  Hence  for  a 
temperature  change  of  20°C,  the  change  in  the  index  of  refrac¬ 
tion  amounts  to  approximately  6  x  10“3.  The  contribution  due  to 
fluctuations  in  the  dopant  concentration  is  estimated8  to  be 
1 0-2  for  a  concentration  of  the  order  of  10 18  cm-3  at  900  nm. 

If  a  diode  laser  array  is  produced  in  a  single  crystal,  then 
cleavage  along  a  crystal  plane  yields  a  typical  error9  in  the 
crystal  length  of  the  order  of  AL  =  3  ym.  If  one  assumes  that 
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the  index  of  refraction  of  the  gain  medium  is  3.6. ,  the  gain 
length  is  300  ym,  and  the  wavelength  is  850  run,  Eq.  (3)  yields 

A<(>  =  1.1  x  10 ^  rad.  (4) 

t 

which  is  a  significant  phase  change  between  the  individual  lasers. 

It  is  appropriate  to  determine  the  effects  of  such  phase  shifts 
on  the  far-field  intensity  distribution  or,  equivalently,  the  dis¬ 
tribution  of  the  focused  output  of  a  phased  laser  array.  In  view 
of  recent  experimental  demonstrations10-12  of  diffraction-coupled 
GaAs  diode  lasers,  we  examine  the  problem  by  assuming  a  set  of 
emitters  having  certain  specific  phase  relations.  These  studies 
lead  to  an  understanding  of  the  experimental  demonstrations. 

The  geometry  under  consideration  is  shown  in  Figure  4.  The 


intensity 

distribution  at 

the  focal  spot  f  is 

given  by1 3 
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in  the  Fresnel  limit,  where  <f(x',y')  is  the  input  field  to  the  lens, 

i.e.,  it  is  the  electric  field  distribution  of  the  emitters.  We 

assume  for  simplicity  that  the  array  is  1-Dimensional,  along  the 

x-axis,  and  the  electric  field  is  an  identical  Gaussian  for  all 

th 

emitters.  Hence  the  electric  field  for  the  n —  emitter  is 


(x-x  )2 
n 

2aQ2  i<j>(xn) 
e 


(6) 


t  h 

where  a  is  the  spot  size  and  x  is  the  position  of  the  n—  laser, 
o  th  ^ 

<J>  (x  )  is  the  phase  of  the  n—  laser.  The  evaluation  of  the  inten¬ 
sity  distribution  using  Eq.  (6)  is  straightforward,  and  the  results 
are  now  presented  for  different  sets  of  phase  shift  parameters. 
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Figure  4.  Geometry  for  the  calculation  of 
the  focused  output  of  an  array 
of  emitters. 


The  first  case  to  be  considered  is  shown  in  Figure  5  where  we 

assume  two  emitters  are  in  phase  with  each  other.  The  results  to 

be  shown  have  the  separation  between  lasers  normalized  to  the  focal 

spot,  i.e.  ,  the  normalization  constant  is  a  =  Af/2Tra  .  a  is 

e  o  o 

taken  to  be  3  ym.  Figure  6  represents  the  focused  intensity  when 
the  separation  between  the  emitters  is  5  ym.  As  the  separation 
increases  to  10  ym,  one  observes  that  additional  side  lobes  appear 
as  shown  in  Figure  7.  This  behavior  arises  from  simple  diffraction 
and  interference  of  the  electromagnetic  waves.  The  existence  of  the 
central  lobe  in  the  output  intensity  is  an  intrinsic  property  of  the 
in-phase  assumption. 
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13701-2 


Figure  5.  Input  field  distribution  for  two 
in-phase  emitters. 
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Figure  6.  Output  intensity  for  the  in-phase  case 
where  the  separation  between  emitters 
is  5  ym.  Transverse  distance  is  meas¬ 
ured  in  units  of  (Xf/2ira0). 


22 


INTENSITY 


13701-4 


TRANSVERSE  DISTANCE 


Figure  7.  Output  intensity  for  the  in-phase  case 
where  the  separation  between  emitter 
is  10  yin.  Transverse  distance  is  meas¬ 
ured  in  units  of  (Af/2Trao). 
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In  contrast  to  the  in-phase  situation.  Figure  8  considers 
the  case  of  two  emitters  which  are  180°  out  of  phase.  The  focused 
output  intensity  distribution  is  shown  in  Figure  9  for  the  case 
where  the  separation  between  the  emitters  is  5  ym.  The  absence  of 
the  central  lobe  (in  comparison  to  the  in-phase  case)  is  due  to  the 
destructive  interference  of  the  two  electromagnetic  waves  which 
are  180°  out  of  phase  with  respect  to  each  other.  Increasing  the 
separation  of  10  ym  yields  additional  side  lobes  due  to  diffraction 
effects . 


Figure  8.  Input  field  distribution  for  two 
180°  out-of-phase  emitters. 

Next,  we  consider  cases  representative  of  the  experiments 
reported  in  the  literature.10-12  We  assume  10  emitters,  each 
having  a  spot  size  of  the  order  of  3  ym.  Figure  10  illustrates 
the  input  field  distribution.  First  we  shall  consider  the  case 
where  all  emitters  are  in-phase.  The  output  intensity  distribu¬ 
tions  are  shown  in  Figures  11  and  12.  Again,  the  central  lobe  is 
very  apparent  in  both  figures  and  is  due  to  the  constructive  nature 
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Figure  9.  Output  intensity  for  the  out-of-phase 
case  where  the  separation  between 
emitters  is  5  ym.  Transverse  distance 
is  measured  in  units  of  (Xf/2Tra0). 
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Figure  10. 


Input  field  distribution  for  10  emitters 
which  are  in-phase  with  each  other. 
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Figure  11.  Output  intensity  for  the  case  when  the 
separation  between  in-phase  emitters 
is  4  ym.  Transverse  distance  is 
measured  in  units  of  (Af/27Ta  )  . 
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Figure  12.  Output  intensity  for  the  case  when  the 
separation  between  in-phase  emitters 
is  10  pm.  Transverse  distance  is 
measured  in  units  of  (A.f/2irao)  . 
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of  the  interference  of  the  eletromagnet ic  waves.  As  the  separation 
between  emitters  is  extended  from  4  pm  (Figure  11)  to  10  ym  (Fig¬ 
ure  12),  the  additional  side  lobes  appear  due  to  (higher  order) 
diffraction.  For  the  purpose  of  comparison,  we  have  plotted  in 
Figure  13  the  output  intensity  profile  for  the  case  when  the 
separation  between  the  emitters  is  10  ym,  but  the  phases  between 
them  are  completely  random.  The  case  where  the  10-emitter  system 
has  180°  phase  shift  between  adjacent  lasers  is  shown  in  Figure  14. 
The  output  intensities  are  shown  in  Figure  15  for  the  case  when  the 
separation  between  emitters  is  4  ym.  Increasing  the  separation  to 
10  ym  gives  rise  to  the  enhancement  of  certain  side  lobes  as  shown 
in  Figure  16.  The  introduction  of  random  phases  yields  a  completely 
random  intensity  distribution  as  shown  in  Figure  17.  It  is 
interesting  to  note  that  the  ten-element  Xerox  experiment14*  was 
indeed  producing  a  diffraction-coupled  phased-array  laser.  The 
output  intensity  of  the  experiment  is  given  in  Figure  18,  which 
shows  good  qualitative  agreement  with  the  simple  analytical 
modeling  represented  by  Figure  16. 


2.2  LOCKING  OF  MODULAR  LASERS  VIA  PHASE  CONJUGATION 

As  was  implied  in  Figure  1,  multiple  independent  gain  modules 
can  be  locked  in  phase  and  frequency  without  mutual  coupling  by 
employing  a  common  four-wave  mixer  at  one  end  of  the  resonator. 
Alternatively,  one  may  substitute  a  set  of  three  four-wave  mixers 
with  a  common  pump  system  for  the  single  mixer.  Each  separate 
resonator  in  this  system  forms  a  stable,  phase-conjugate  resonator 
in  which  intra-resonator  aberrations  are  internally  compensated 
if  the  output  mirrors  are  aberration  free.  Further,  if  the  output 
mirrors  are  precisely  aligned  with  each  other,  i.e.,  if  6L^  =  SL^ 

=  6L^  in  Figure  19,  then  the  output  fields  are  also  in  phase. 
Although  the  illustrated  output  mirrors  are  flat  segments,  one  could 
substitute  curved  mirrors,  as  appropriate  to  a  confocal  resonator, 
for  example,  with  the  same  basic  locking  performance. 
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Figure  13 .  Output  intensity  of  the  10-emitter  system 

when  the  phases  between  emitters  are 

completely  random.  The  separation  between 

emitters  is  10  pm.  Transverse  distance  is 

measured  in  units  of  (Af /2?ra  )  . 
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Figure  14.  Field  distribution  for  the  case  of 
180°  phase  shift. 
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Figure  15.  Intensity  distribution  for  the  case  when 
the  adjacent  emitters  are  180°  out-of¬ 
phase  and  the  separation  between  emitters 
is  4  ym.  Transverse  distance  is  measured 
in  units  of  (Af/27Tao)  . 
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Figure  16.  Intensity  distribution  for  the  case  when  the 
adjacent  emitters  are  180°  out-of-phase  and 
the  separation  between  emitters  is  10  urn. 
Transverse  distance  is  measured  in  units  of 
( Af /2TTaQ)  . 
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Figure  17.  Intensity  distribution  for  the  case  when 
the  10-emitter  system  has  random  phases 
among  all  emitters.  Transverse  distance 
is  measured  in  units  of  (Xf/2iTao)  . 
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ANGLE,  DEGREES 

Figure  18.  The  output  intensity  of  the  Xerox 
phased  array  diode  laser.  (From 
Reference  14).  Curves  corresponding 
to  different  output  powers  have  been 
separated  vertically  for  clarity. 
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Figure  19.  Geometry  of  the  phase  locking  analysis 

with  a  four-wave-mixer  conjugation  lock¬ 
ing  technique. 
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In  order  to  analyze  the  requirements  needed  to  achieve  a  high 
quality  output  beam,  consider  the  field  at  an  arbitrary  plane, 
located  at  position  z',  inside  of  each  resonator,  as  illustrated 
in  Figure  19.  The  three  output  coupling  mirrors  are  assumed  to  be 
spaced  from  this  plane  by  variable  distance  <SLn-  Ignoring  trans¬ 
verse  wavefront  variations,  the  field  incident  on  plane  z'  from 
the  right  has  a  complex  amplitude  where 

-i<t>  (z1  ) 

A'  =  Ae  n  ,  (7) 

n 

where  A  is  a  scaler  and  <j>  is  a  phase  shift  which  is  to  be  deter¬ 
mined.  After  reflecting  off  the  output  coupling  mirror  and  return¬ 
ing  to  plane  z' ,  the  field  amplitude  is: 


An  =  A  \/^exp  ("i(<f)n  (Z,)  "  2k6Ln  +  VH  '  (8) 


where  <J>  is  the  mirror-induced  phase  shift  and  is  the  power 
reflectivity  of  the  output  mirror. 

After  amplification  and  reflection  from  the  conjugator, 
followed  by  reamplification,  the  returning  field  at  plane  z’  has 
a  complex  amplitude  A' 1 ' ,  where 


A' ' '  =  A  JR  R  G 
n  v  m  c  s 


exp 


i ($n (z 1 )  - 


2k6L  +  <f>  +  d>  . 

n  Ym  Ynl 


4>n2> 


(9) 


where  G  is  the  steady  state  round-trip  power  gain  of  the  amplifier 
s 

Rc  is  the  conjugator  reflectivity,  and  <J>  ^  and  §n2  are  *-he  P^ases 
of  pumps  1  and  2,  respectively.  After  many  iterations,  the  oscilla 
tion  condition  requires  that  A' ' '  =  A^,  so  that 

Jr  R  G  =  1  ,  (10) 

v  c  m  s 
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and 


<f>n(z')  =  1/2  [-2k6Ln  +  4>m  +  <Pnl  +  <pn2]  .  (11) 


Thus,  if  the  output  mirrors  are  aligned  so  that  all  6Ln  are  equal, 
and  if  c|>  is  the  same  for  all  mirrors,  and  if  the  mixer  pumping  is 
coordinated  so  that  <}>n^  +  <p^2  is  the  same  for  all  four-wave  mixers, 
then  all  three  phases  <j>  (z1.)  are  equal.  That  is,  the  output  wave- 

front  is  coplanar.  The  fact  that  the  output  phase  front  is  not 
flat  if  the  end  mirrors  are  offset  from  each  other  is  analogous  to 
the  situation  on  a  standard  PCR  with  an  aberrated  out-coupling 
mirror.  In  this  case,  it  has  been  shown  that  the  phase  of  the  out- 
coupled  wave  matches  the  distortion  of  the  mirror.15  The  impact 
of  this  analysis  on  coherent  phasing  of  multiple  GaAs  diodes  is 
either  that  one  well-figured  mirror  must  be  used  for  the  whole 
array  or  that  the  separate  end  mirrors  must  be  precisely  aligned 
to  much  better  than  a  wavelength. 

Preliminary  three-dimensional  computer  calculations  of  a 
multiple  diode  PCR  show  that  the  output  of  such  a  resonator  is 
virtually  independent  of  refractive  index  and  length  variations 
of  the  individual  diodes.  Figure  20  shows  the  amplitude  and  phase 
of  a  converged  solution  to  the  resonator  shown  in  Figure  19  at 
various  locations  in  the  resonator.  Parameters  used  for  this 
sample  calculation  are: 


Wavelength 
Diode  cross  section 
Gain-length  product 
Saturation  intensity 
Distance  to  conjugator 
Output  coupling 


840  nm 

13.50  x  40  ym 
0.50 

40.00  kW/cm^ 
500.00  ym 
50  percent 
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(b)  INTRA-CAVITY  INTENSITY  AND  PHASE 


(c)  OUTPUT  FAR-FIELD  INTENSITY 

Figure  20.  Computer  calculations  showing  intensity 
and  phase  at  different  locations  for  a 
three-diode  phase  conjugate  resonator. 
Phase  distortions  arise  from  different 
diode  lengths. 
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Despite  diode-to-diode  phase  variations  of  more  than  half  a  wave, 
the  output  is  flat.  For  a  single  diode  of  the  same  characteristics, 
the  output  amplitude  and  phase  are  given  in  Figure  21.  The  total 
power  extracted  from  one  diode  is  17  mW  while  from  three  diodes 
47  mW  is  extracted.  The  loss  in  extraction  efficiency  per  diode 
in  the  array  is  less  than  8  percent.  For  multimode  waveguide 
diodes  the  loss  should  be  even  smaller. 

Calculations  have  also  been  performed  with  a  conventional 
mirror  in  place  of  the  PCM.  Even  with  no  phase  variations  between 
the  diodes,  the  output  is  poor  as  shown  in  Figure  22  which  gives 
the  far-field  irradiance  pattern  of  this  calculation  after  twenty 
iterations.  The  pattern  is  clearly  inferior,  and  furthermore  the 
power  is  less  than  1  mW  at  this  iteration. 

2.3  BEAM  COMPACTING  SYSTEMS  FOR  HIGH  QUALITY  FAR-FIELD 

PERFORMANCE 

Arrays  of  laser  diodes  must  necessarily  be  packed  with 
appreciable  spacing  between  the  diodes  (relative  to  the  diode 
size)  in  at  least  one  dimension  -  the  narrow  dimension  of  the 
diode.  In  order  to  cool  the  diodes  better,  it  may  also  be 
important  to  space  the  diodes  in  the  direction  of  the  long 
dimension  of  the  diode.  Since  far-field  power  combining 
appears  to  be  the  most  expedient  and  efficient  technique  for 
combining  the  power  outputs  of  N  diodes,  it  is  important  to 
examine  the  impact  of  a  widely  spaced  set  of  radiation  sources 
on  the  far-field  distribution.  This  is  an  extremely  important 
consideration  in  diode  arrays  intended  to  drive  nonlinear  opti¬ 
cal  devices,  since  such  systems  typically  focus  the  available 
power  into  the  nonlinear  crystal  in  order  to  generate  sufficient 
optical  intensity  to  achieve  efficient  doubling.  The  field 
distribution  at  the  focus  is  essentially  a  far-field  pattern. 
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Figure  21. 


Output  intensity  and  phase  of  a  single 
diode.  (Computer  simulation) 
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Figure  22. 


Far-field  intensity  of  a  3-diode  laser 
with  a  conventional  feedback  mirror. 


Assuming  that  the  diode  array  is  arranged  on  a  periodic 
lattice,  as  viewed  end-on,  with  appreciable  spacing  between  ele¬ 
ments,  then  the  radiating  source  resembles  a  diffraction  grating. 

It  is  well  known  that  the  far-field  pattern  of  such  a  uniphase 
source  also  has  a  quasi-periodic  character  with  a  period  determined 
by  the  inverse  of  the  spacing  between  diodes,  as  exhibited  in 
Figure  23b,  which  is  taken  from  Reference  16.  This  pattern  con¬ 
sists  of  a  major  lobe  on  axis  plus  a  large  number  of  similarly 
shaped  secondary  or  "grating"  lobes  with  the  nearby  lobes  of 
comparable  size  to  the  main  lobe  (for  large  percentage  diode 
spacings) .  Unfortunately,  for  the  diode  spacings  required  in  our 
modular  laser  application,  most  of  the  far-field  power  is  contained 
in  the  composite  of  the  grating  lobes  rather  than  in  the  main  lobe. 
In  the  practical  applications,  the  power  in  these  grating  lobes 
will  effectively  be  lost. 
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Figure  23.  Near-  and  far-field  patterns  of  well  aligned 
and  poorly  aligned  arrays  (after  Ref.  16). 
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Randomizing  the  position  of  the  diodes  destroys  the  periodicity 
of  the  source,  and  of  the  far-field  pattern  as  well,  such  that  the 
side  lobe  peaks  are  less  intense  than  with  the  periodic  array.  This 
is  dramatically  illustrated  in  Figure  23d.  Unfortunately,  the  total 
power  which  is  lost  from  the  main  lobe,  given  the  same  average  diode 
spacing,  is  unchanged  by  the  randomization.  This  loss  is  determined 
almost  entirely  by  the  area  thinning  factor  of  the  diode  array, 
which  is  defined  as  the  ratio  of  the  total  diode  effective  radiating 
area  to  the  area  of  the  box  encompassing  the  diode  (end  face)  array. 
The  effective  intensity  is  directly  proportional  to  the  thinning 
factor.  Thus  for  a  thinning  factor  of  50%,  half  of  the  radiated 
power  is  contained  in  the  side  lobe  structure  and  is  thereby  lost 
to  the  main  lobe,  insofar  as  the  diodes  maintain  their  average 
spacing  under  randomization. 

The  only  workable  solution  to  this  problem  for  a  given  diode 
spacing  is  to  use  mirrors  (Figure  24a)  or  lenses  (Figure  24b) , 
or  some  combination  of  these  techniques,  to  redistribute  the  diode 
near-field  pattern  in  such  a  way  that  a  more  compact  radiating 
source  is  produced.  The  mirror  systems  illustrated  in  Figure  24a 
tend  to  produce  unequal  path  differences  which  could,  in  principle, 
be  compensated  with  optical  glass  slabs  of  differing  length.  How¬ 
ever,  precision  brute  force  path  matching  seems  extremely  difficult. 
Similar  path  matching  problems  occur  for  the  lens  arrays  shown 
in  Figure  24b.  One  major  advantage  of  using  a  phase  conjugate 
resonator  (PCR)  in  this  diode  coupling  application  is  that  these 
path  length  differences  are  automatically  cancelled  out  if  the 
beam  combining  optics  are  placed  within  the  PCR. 

2 . 4  NORMAL  MODES  OF  A  MULTIPLE-ELEMENT  PLANAR  WAVEGUIDE 

This  section  presents  calculations  of  the  normal  modes 
of  a  set  of  coupled  emitters.  This  analysis  is  relevant  to 
coupling  schemes  such  as  the  use  of  a  phase  conjugate  resona¬ 
tor  or  evanescent-wave  coupling.  In  such  schemes,  the  output 
of  a  phased-array  laser  system  can  be  considered  in  terms  of 
the  superposition  of  the  normal  modes  of  their  respective 
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Figure  24.  Beam  compacting  systems  for  diode-laser 

arrays.  The  feedback  mirrors  may  be  moved 
to  the  left  of  the  beam-compacting  elements 
in  order  to  compensate  for  optical  path 
differences . 
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waveguide  system.  This  picture  has  certain  limitations.  It 
cannot  be  rigorously  applied  to  a  saturated  gain  medium  due  to  the 
nonlinear  nature  of  the  index  of  refraction;  i.e.,  the  detuning 
of  the  stimulated  photon  from  the  transition  resonance  gives  rise  to 
an  effective  nonlinear  susceptibility.  The  saturated  regime  poses 
an  extremely  difficult  mathematical  problem.  Only  recently, 
numerical  solutions  have  been  obtained  that  provide  a  good  physical 
insight  into  the  role  of  the  saturated  index  of  refraction  on  the 
transverse  mode  of  a  single  gain  medium. 1 7  However  no  work  has 
been  done  on  an  array  of  laser  oscillators.  Hence,  we  would  like 
to  provide  a  first  step  toward  the  understanding  of  the  behavior 
of  phased-array  lasers  by  analyzing  the  multiple  waveguide  problem. 
This  approach  is  valid  for  the  linear  and  homogenous  gain  regime. 

By  homogenous  gain  we  shall  imply  that  the  gain  coefficient  is 
constant  for  each  laser  and  is  the  same  for  all  lasers. 

The  problem  under  consideration  is  illustrated  in  Figure  25. 

The  starting  point  of  the  analysis  is  the  scalar  wave  equation 


0 


(12) 


in  the  harmonic  approximation.  The  scalar  field  can  be  written 
as 


E 


(13) 


such  that  the  slowly  varying  envelope  approximation 


(14) 


is  satisfied.  Equation  (14)  implies  that  both  the  amplitude  and 
phase  of  the  complex  envelope  $  vary  slowly  compared  to  the 
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Figure  25.  Geometry  for  the  multiple  planar 
waveguide.  ng  and  n^  are  the 
indices  of  refraction. 


effective  wavelength  3  of  the  waveguide.  Application  of 
Eqs.  (13)  and  (14)  to  the  wave  Eq.  (12)  yields 


2i3 


M 

3z 


+ 


0 


(15) 


which  is  just  the  quasi-optical  or  paraxial  wave  equation.  Eq.  (15) 
has  been  the  starting  point  of  numerous  numerical  simulations  of 
propagation  in  optical  fiber  systems.18  The  normal  modes  u  of 
the  system  are  defined  to  be  the  solutions  of 


(16) 
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Any  general  wave  can  be  written  as  a  superposition  of  the  normal 
modes 


*(x,y)  =  Aa  (z)  ua  (x)  (17) 

a 


The  coeffieicnt  A^fz)  can  be  computed  in  terms  of  the  initial 
field  distribution.  Hence  the  essence  of  this  section  is  to  calcu¬ 
late  the  normal  modes  of  the  waveguide  as  a  function  of  the  trans¬ 
verse  dimension  of  each  guide.  Since  Eq.  (16)  is  a  second-order 
differential  equation,  we  note  that  the  condition 

2 

^2  n2  -  e2  >  0  (18) 

c 

is  required  for  the  existence  of  confined  modes,  and 


32  <  0 


(19) 


is  required  for  the  existence  of  unconfined- modes .  Since  the  indices 
of  refraction  are  chosen  such  that  for  n.^  the  situation  correspond¬ 
ing  to  Eq.  (18)  is  satisfied,  while  for  nQ  condition  Eq.  (19)  is 
satisfied,  we  are  dealing  with  a  configuration  shown  in  Figure  26. 
Boundary  conditions  must  then  be  satisfied;  they  are  given  by 


u  (x) 

continuous  at 
the  boundary 

du 

dx 


(20) 
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Figure  26.  Relative  size  of  the  indices  of 
refraction  for  the  distinct 
waveguides . 


In  the  same  spirit  as  in  Section  2.1,  we  shall  build  our 
knowledge  by  solving  relatively  simple  problems  first.  In  this 
context,  we  consider  the  normal  modes  of  a  single  planar 
waveguide  whose  geometry  is  shown  in  Figure  27.  We  shall  assume 
that  the  distribution  of  the  magnitudes  of  the  index  of  refrac¬ 
tion  are  such  that  n^  >  n^ ,  as  shown  in  Figure  28.  Hence,  the 
effective  wavenumbers,  are  defined  such  that 


k 


2 

1 


e2  >  o. 


I X I  <  L 


$2  <  0  , 


I  x  |  >  L 


(21) 
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Figure  27.  Geometry  of  the  single 
planar  waveguide. 
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Figure  28.  Relative  sizes  of  the  index  of 

refraction  for  the  single  planar 
waveguide  case. 
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The  solutions  can  be  obtained  by  inspection.  They  are 


u  (x) 

u  (x) 


u  (x) 


k_.x  ^  _ 

=  A  e  0  ,  x  <  -L 

{B  cos  k,x 

1  l*|  <  L 

C  sin  k^x 

=  D  e"k0X,  x  >  L 


(22) 


Applying  the  boundary  conditions  at  x  =  L,  one  finds  the  following 
dispersion  relations 


,  T  0  even 

tan  k.,L  =  ^ — ,  -  . 

1  k^  solution 

.  __  i_  T  _  kl  odd 
tan  kjL  solution 


(23) 


By  even  solution,  we  mean  the  cosine  function  solutions  in  the 
region  |x|  <  L,  while  odd  solutions  are  represented  by  sine  func¬ 
tions  in  the  same  region.  The  dispersion  relations  determine 
precisely  the  number  of  modes  available  for  a  specific  transverse 
dimension  L.  Figures  29  and  30  show  the  plot  of  the  lowest  order 
for  the  even  and  odd  modes,  respectively. 

Next,  we  consider  the  normal  modes  of  the  double  planar 
waveguide  whose  gemoetry  is  shown  in  Figure  31.  The  relative 
indices  of  refraction  are  shown  in  Figure  32.  Following  the  same 
procedure  as  outlined  above  the  solutions  are 
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Figure  29.  Lowest  order  even  mode. 
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Figure  31.  Geometry  for  the  double 
planar  waveguide. 
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Figure  32.  Relative  size  of  the  index  of  refraction 
as  a  function  of  the  transverse  dimension 
of  the  waveguide  system. 
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The  following  dispersion  relations  are  found  in  the  same  manner 
as  before. 


a)  Even  symmetric  solutions,  or  odd  antisymmetric  solutions 


tan 


=tanh  kgL 


(25) 


tan  k^  (£  +  L) 


b)  Even  antisymmetric  solutions,  or  odd  symmetric  solution. 


tan  k^£ 


-coth  kg& 


(26) 


tan  k^  (£  +  L) 


Figure  33  shows  a  plot  of  the  dispersion  relation  as  a  function 
of  the  product  of  3  and  the  waveguide  dimension  L  for  the  case  of 
even-symmetric  or  odd-antisymmetric  solutions.  Figure  33  shows 
two  curves  which  represent  the  right  and  left  hand  sides  of 
Eq.  (25) .  The  intersection  points  are  the  allowed  values  of  3 
for  fixed  lengths  SL  and  L.  For  example,  if  3L  =  67.6,  one  can 
count  that  there  exist  five  allowed  normal  modes. 

The  next  case  to  be  considered  is  the  triple  planar  waveguide 
system  illustrated  in  Figures  34  and  35.  Again  the  solutions  can 
be  written  in  a  similar  manner  as  before,  and  the  application  of 
boundary  conditions  yields  the  dispersion  relations.  For  simplicity 
we  shall  consider  the  even-symmetric  solutions.  The  dispersion 
relation  is  given  by 


k 


0 


tan  k2  (21  +  3^f)  tanh  kQ 


(21  +  |) 


(27) 
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1 


tank^  (21  +  If) 


tanh 
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BETA-WAVEGUIDE  LENGTH 

Figure  33.  Dispersion  relation  for  the  even -symmetric 
or  odd-antisymmetric  solutions  of  the 
double  planar  waveguide  system. 
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Figure  34.  Geometry  for  the 

triple  planar  wave¬ 
guide  system. 
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Figure  35.  Relative  magnitude  of  the  indices  of 
refraction  as  a  function  of  the 
transverse  dimension  of  the  wave¬ 
guide  system. 


Equation  (27)  is  plotted  in  Figure  36  showing  two  curves  corres¬ 
ponding  to  the  right-  and  left-hand  sides  of  the  equation.  Iden¬ 
tical  interpretation  can  be  assigned  to  the  significance  of  the 
crossing  or  intersection  points.  The  dispersion  for  the  even- 
antisymmetric  mode  is  given  by 


tan  k1  {21  +  3j) 
tan  k^  {21  +  Ijf) 


coth  k 


0 


{21  +  |) 


coth 


koL 


tan 


k^L 


(28) 


and  a  similar  curve  can  be  obtained  from  Eq.  (28) . 

The  double  and  triple  planar  waveguide  systems  form  the 
basis  for  obtaining,  by  inspection,  the  mode  structure  of  the  2N 
and  2N+1  planar  waveguide  systems,  respectively.  The  conclusion 
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Figure  36.  Dispersion  relation  for  the  case  of 
even  symmetric  solutions  in  the 
triple  planar  waveguide  system. 
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that  we  can  draw  is  that  there  are  intrinsically  4  types  of 
normal  modes;  even  symmetric  and  antisymmetric,  and  odd  symmetric 
and  antisymmetric.  Any  input  field  can  then  be  expressed  in  terms 
of  these  modes.  The  analysis  of  the  phased-array  diode  laser 
system  can  benefit  from  the  knowledge  that  the  output  field  is  no 
more  than  the  superposition  of  these  modes  with  the  appropriate 
coefficients.  The  coefficients  are  obtained  from  the  gain  equations 
for  the  laser  oscillation. 

The  analysis  of  this  section  has  been  used  to  describe  the 
nature  of  the  normal  modes  that  will  arise  in  a  coupled  multi¬ 
element  diode  laser  array.  This  analysis  is  relevant  for  many 
coupling  schemes,  including  the  evanescent-wave  approach  repre¬ 
sented  by  References  10-12,  and  the  phase-conjugate  resonator 
approach  which  is  the  subject  of  the  present  program.  These 
analytical  results  will  ultimately  be  helpful  in  understanding 
which  modes  should  be  selected  by  some  means  in  order  to  achieve 
the  optimum  far-field  intensity  distribution. 
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SECTION  3 


CANDIDATE  PHASE  CONJUGATE  MIRRORS 


The  use  of  a  phase  conjugate  resonator  to  couple  multiple 
diode  gain  elements  was  introduced  in  Section  1.  In  principle, 
aberrations  arising  from  differences  in  the  individual  gain  ele¬ 
ments  are  automatically  compensated  by  the  phase  conjugate 
mirror  (PCM) .  This  section  considers  what  PCM  concepts  are 
compatible  with  the  wavelength  range  and  power  levels  of  the 
multi-element  GaAlAs  diode  laser.  Specifically,  we  have  con¬ 
sidered  phase  conjugation  by  the  process  of  stimulated  Brillouin 
scattering  (SBS) .  This  approach  avoids  the  requirement  for 
external,  high  quality  pump  lasers  necessary  to  construct  a  PCM 
based  on  four-wave  mixing.  However,  because  of  the  <  1  GHz  fre¬ 
quency  shift  imparted  to  the  laser  radiation  as  part  of  the  SBS 
process,  this  approach  can  only  be  utilized  in  pulsed  applications 
where  the  bandwidth  of  the  gain  medium  is  broad  enough  to  accomo¬ 
date  the  accumulated  frequency  shift.  For  a  4-nm  bandwidth,  a 
1 0-cm-long  resonator  could  operate  at  pulse  lengths  up  to  1  ysec. 

Stimulated  Brillouin  scattering  is  a  threshold  process: 
only  if  the  incident  light  intensity  exceeds  a  minimum  value 
will  appreciable  reflection  occur.  Hence,  a  crucial  condition 
is  the  threshold  intensity  for  SBS  to  arise.  The  objective 
of  this  section  is  to  calculate  such  a  threshold  intensity  for 
two  types  of  optical  waveguides.  One  is  filled  with  carbon 
disulfide  (CS^)  and  the  other  one  is  composed  of  fused  quartz.  The 
expression  for  the  threshold  intensity  is  given  by19 
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where  c  is  the  velocity  of  light,  T  is  the  bulk  modules,  e  is  the 
permeability  of  free  space,  n1  and  n2  are  the  indices  of  refractior 
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at  the  pump  and  SBS  frequencies,  respectively,  y  is  the 
electrostriction  coefficient,  is  the  wavevector  of  the  pump; 
k-2  is  the  wavevector  of  the  acoustic  wave,  is  the  attenuation 
length  of  the  pump  and  Lg  is  the  attenuation  length  of  the  acoustic 
wave.  For  the  case  of  CS2,  one  finds  that  the  threshold  intensity 
is  approximately  given  by 

t  _  a  .7 r“  i  a5  v; atts  , _  A . 

Ith  =  0-75  x  10  - 2~  (30) 

cm 

The  threshold  energy  for  SBS  can  then  be  calculated  as 

E.^  =  0.75  x  10^  watts  x  Area  x  T  (31) 

cm 

where  x  is  the  pulse  duration,  and  the  area  is  the  input  cross- 
section  at  the  entrance  of  the  multimode  waveguide  system.  Fig¬ 
ure  37  shows  a  plot  of  the  threshold  energy  versus  optical  fiber 
diameter  for  distinct  values  of  x.  The  optical  fiber  is  assumed 
to  be  filled  with  liquid  CS2-  The  wavelength  chosen  for  this 
calculation  is  700  nm.  However,  CS2  has  a  large  bandwidth  for 
SBS  and  should  be  able  to  accommodate  all  GaAlAs  diode  laser 
wavelengths.  In  the  case  of  Si02  optical  fiber  one  finds  that 
the  threshold  power  is20 
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(32) 


where  A  is  the  cross  section  of  the  optical  fiber,  g  is  the  SBS 
gain  and  is  the  damping  length  of  the  electromagnetic  wave. 

Experimental  measurements20  yield  g  =  4.3  x  10-9  cm/Watt  and 
Lef£  =  2  x  10  cm.  The  threshold  energy  is  then 

Eth  =  1.72  x  10^  x  Area  x  x  (33) 
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Figure  37.  SBS  threshold  energy  versus  optical  fiber 
diameter  for  several  pulse  lengths.  The 
optical  fiber  is  assumed  to  be  filled  with 
liquid  CS2- 
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Figure  38  shows  a  plot  of  Eth  vs  fiber  diameter  for  the  Si02 
fiber.  These  results  indicate  that  the  threshold  energy  require¬ 
ment  for  the  generation  of  phase  conjugate  waves  by  SBS  is  rather 
low.  The  smallness  of  the  value  is  due  to  two  factors;  long 
attenuation  length  for  the  electromagnetic  wave,  and  small  input 
cross-section  for  the  optical  fibers. 


Figure  38.  SBS  threshold  energy  versus  optical  fiber 
diameter  for  several  pulse  lengths.  The 
optical  fiber  is  made  of  SiC>2. 
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SECTION  4 


HIGH  EFFICIENCY  FREQUENCY  DOUBLING  CONCEPTS 

4 . 1  INTRODUCTION 

In  this  section  we  show  that,  if  concepts  for  the  co-phasing 
of  diode  gain  media  can  be  successfully  implemented,  then  fre¬ 
quency  doubling  with  efficiencies  consistent  with  the  requirements 
of  practical  military  applications  appear  possible.  This  calcu¬ 
lation  is  based  on  published  results21  describing  the  successful 
(but  low  efficiency  <0.2  percent)  demonstration  of  frequency 
doubling  using  the  crystal  KNbO^,  driven  by  a  GaAlAs  diode  laser 
operating  at  860  nm.  By  scaling  these  results  and  using  the  high 
peak  power  and  narrow  linewidth  properties  anticipated  from  the 
proposed  diode  phased  arrays,  it  appears  that  as  few  as  thirteen 
commercial  diode  lasers  operating  coherently  will  provide  suffi¬ 
cient  peak  power  to  achieve  40  percent  doubling  efficiency.  In 
addition,  the  following  analysis  indicates  that  high  efficiency 
will  be  achievable  at  kilowatt  power  levels,  providing  the 
doubling  process  can  be  scaled  without  sacrificing  efficiency. 
Hughes  feels  this  can  be  accomplished  with  a  minimum  number  of 
nonlinear  crystals  based  on  considerable  experimental  and  theo¬ 
retical  research  on  the  problems  involved  in  scaling  frequency 
doubling  devices  to  high  average  power. 

In  this  section  we  first  lay  the  groundwork  for  predicting 
efficiency  using  well  established  principles.  These  principles 
are  then  used  to  select  a  doubling  crystal,  where  we  find  that 
KNb03  is  the  best  candidate  for  the  present  application,  although 
other  materials  are  also  identified.  We  then  describe  the  scaling 
concepts  which  lead  to  the  high  efficiency  (40  percent)  doubling 
potential.  A  discussion  of  the  solutions  to  frequency-doubler 
scaling  problems  is  then  presented. 
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4.2  NONLINEAR  CRYSTAL  SELECTION 


This  section  reviews  the  physical  principles  that  ensure 
efficient  frequency  doubling  in  the  absence  of  the  severe  thermal 
effects  that  arise  in  scaling  to  high  average  power.  (Power 
scaling  is  discussed  in  Section  4.4).  Consideration  of  these 
principles  leads  to  a  set  of  selection  criteria  for  a  nonlinear 
crystal  appropriate  for  the  present  application.  Several  crystal 
candidates  are  described.  In  the  following  section  one  material 
is  selected  and  baseline  performance  is  predicted  using  this 
material . 

Frequency  doubling  has  been  studied  for  twenty  years;  the 
basic  analytical  equations  describing  the  process  are  readily 
available,22  and  show  that  the  doubling  efficiency  for  a  funda¬ 
mental  wave  at  a  frequency  w  with  its  power  P1  uniformly  distri¬ 
buted  over  an  area  A  is  given  by 


3/2  2  , 2  2  , 

i  to  d  L  | 

("sin  (AkL/2)"| 

1  3 

n  1 

L  AkL/2  J 

where  d  is  the  appropriate  nonlinear  coefficient  for  a  crystal 
with  an  interaction  length  L  and  refractive  index  n,  and 

Ak  =  k2  “  2kx  (35) 

is  the  difference  between  the  second- harmonic  wavevector  and 
twice  the  fundamental  wavevector  k-.  Equation  (34) ,  which 
neglects  depletion  of  the  fundamental  beam,  shows  that  doubling 
efficiency  is  maximized  by  selecting  a  material  with  a  high  non¬ 
linear  coefficient,  a  long  interaction  length,  and  by  arranging 
for  phase  matching  to  occur  (Ak  =  0) .  Phase  matching  occurs 
when  the  fundamental  and  second- harmonic  beams  propagate  through 
the  crystal  with  the  same  phase  velocity. 
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The  technique  to  achieve  phase  matching  utilizes  the  natural 
birefringence  of  anisotropic  crystals.  With  a  suitable  choice  of 
the  polarization  and  the  propagation  direction  of  the  funda¬ 
mental  wave,  there  often  exists  a  different  polarization  for  the 
second  harmonic  that  allows  propagation  at  the  same  phase  velocity 
as  the  fundamental.  When  considering  practical  applications  in 
which  high  doubling  efficiency  is  imperative,  however,  this 
"critical"  phase  matching  is  undesirable  for  two  reasons.  First, 
as  the  extraordinary  wave  propagates  through  the  crystal,  its 
direction  of  power  flow  differs  from  the  direction  of  its  wave- 
vector.23  This  leads  to  a  gradual  walk-off  of  the  extraordinary 
wave  from  the  ordinary  wave,  and  severely  limits  the  effective 
interaction  length  of  the  nonlinear  crystal. 2 4  Second,  the 
required  angular  tolerance  on  the  propagation  direction  of  the 
fundamental  can  typically  be  as  small  as  1  mrad  for  a  1  cm  inter¬ 
action  length,  thereby  limiting  the  allowed  beam  divergence  of 
the  fundamental  wave. 

For  these  two  reasons,  it  is  preferred  to  utilize  crystals 
that  allow  non-critical ,  90  degree  phase  matching.25  In  this 
configuration,  the  beams  propagate  along  a  crystal  axis.  Con¬ 
sequently  the  walk-off  does  not  occur,  and  the  angular  tolerance 
is  greater  by  typically  an  order  of  magnitude.  Often,  noncritical 
phase  matching  may  be  achieved  by  adjusting  the  crystal  tempera¬ 
ture  until  the  refractive  indices  reach  appropriate  values. 
Unfortunately,  this  is  not  always  possible  if  the  refractive 
indices  are  not  sufficiently  temperature  dependent,  or  if  the 
necessary  phase  matching  temperature  lies  outside  practical 
limits.  For  example,  temperatures  <  0°C  are  inconvenient  because 
of  potential  window  frost  problems,  or  the  required  temperature 
might  be  near  a  phase  transition  or  crystal  degradation 
temperature . 

Additional  factors  affecting  crystal  selection  include 
optical  quality  and  uniformity,  transparency  at  both  the  funda¬ 
mental  and  second  harmonic  wavelengths,  resistance  to  surface. 
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bulk,  and  index  damage,  availability  and  damage  resistance  of 
antireflection  coatings,  and  various  physical  properties  such  as 
crystal  stability  and  hygroscopicity . 

Such  considerations  lead  to  the  selection  of  KNbO^  as  the 
optimum  crystal  for  doubling  GaAlAs  radiation.21  Promising 
alternatives  include  KTiOPC>4  (KTP )  ,  recently  the  subject  of 
significant  developmental  effort  supported  by  the  Air  Force, 

BaoNaNbc0..  c ,  and  LiNbO_. 

2  d  lo  3 

4.3  FREQUENCY  DOUBLER  DESIGN 

In  the  work  described  in  Reference  21,  radiation  from  a  diode 
laser  having  0.196  W  of  peak  power  lying  within  the  0.09  nm  phase 
matching  bandwidth  was  focused  into  the  KNbO^  crystal.  The 
observed  second  harmonic  power  of  0.34  mW  (doubling  efficiency 
of  0.18  percent)  is  slightly  more  than  that  expected  from  Eq.  (34). 
In  an  earlier  article,26  the  same  author  reported  a  conversion 
efficiency  of  15  percent  using  a  higher  power  parametric  oscillator 
source  at  the  same  860  nm  wavelength,  but  having  a  narrower  line- 
width  of  0.045  nm.  Those  experimental  results  are  shown  in  Fig¬ 
ure  39.  Once  again,  the  result  was  in  good  agreement  with  the 
predictions  of  Eq.  (34).  It  is  important  to  note  that  no  satura¬ 
tion  effects  for  the  doubling  efficiency  were  evident  for  funda¬ 
mental  intensities  up  to  0.5  MW/cm2,  the  maximum  available  from 
that  pump  laser,  indicating  that  greater  efficiency  is  expected 
for  higher  intensities. 

The  implications  of  this  previous  work  on  the  present  study 
are  significant  in  that  they  show  that  KNb03  can  be  non-critically 
phase  matched  to  produce  second-harmonic  wavelengths  between 
430  nm  and  470  nm,  with  longer  wavelengths  possible  by  raising 
the  crystal  temperature.  In  addition,  the  crystal  behavior  was 
close  to  that  expected  from  theoretical  calculations.  It  is 
appropriate,  therefore,  to  use  those  results  in  estimating  per¬ 
formance  of  such  a  frequency  doubler  when  used  with  a  high  power 
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SECOND  HARMONIC  INTENSITY 


H  562-13 


Figure  39.  Second  harmonic  power  density  as  a  function  of 
fundamental  power  density  for  non- critically 
phase  matched  second  harmonic  generation  at 
room  temperature  in  KNb03-  Crystal  length  is 
5.74  mm.  (From  Reference  26). 
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diode  laser  source  as  an  indicator  of  what  is  expected  when  scaling 
to  large  phased  arrays  of  diodes. 

The  experimental  configuration  of  Reference  21  incorporated 
optical  elements  to  focus  the  radiation  of  a  single  GaAlAs  laser 
into  the  KNbO^  doubling  crystal.  The  optical  design  was  selcted 
to  ensure  that  all  the  radiation  was  just  within  the  acceptance 
angle  of  the  nonlinear  crystal.  In  scaling  up  to  the  number  of 
^iode  elements ,  and  consequently  the  beam  dimensions,  to  increase 
the  power  and  doubling  efficiency,  the  focal  length  of  the  lens 
used  to  focus  the  radiation  into  the  nonlinear  crystal  must 
increase  proportionately,  i.e.,  the  f-number  must  remain  con¬ 
stant,  in  order  not  to  exceed  the  acceptance  angle.  Under  this 
condition  it  can  be  shown  that  if  the  diode  radiation  distri¬ 
butions  are  effectively  compacted  in  such  a  way  as  to  minimize 
"dead"  areas  (see  Section  2.3),  the  intensity  of  the  central 
lobe  of  the  focused  beam  increases  linearly  with  the  number  of 
elements.  This  is  most  easily  understood  by  considering 
the  simple  case  of  linear  array  of  N  plane  wave  sources  of 
dimension  s  separated  by  a  distance  d.  Diffraction  theory27 
predicts  that  in  the  focal  plane  of  a  lens  with  a  focal  length  f 
the  intensity  at  a  distance  x  off  the  optical  axis  is  given  by 


(36) 


where  P  is  the  power  radiated  by  a  single  source,  X  is  the  wave 
length,  and  k  =  2tt/A.  Consider  first  the  case  where  N  =  1;  the 
on-axis  intensity  is 


(37) 
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and  is  determined  by  the  acceptance  angle  of  the  nonlinear 
crystal.  The  on-axis  intensity  for  the  case  of  N  sources  adding 
coherently  is 


V0)  =  y  1,2  (381 

where  f^  is  the  array  focal  length.  In  order  to  maintain  constant 
f-number,  one  requires  that 


(39) 


Equation  (38)  ,  then  becomes 

V°>  “  N(t)(^)=  N  t  V0’  (40) 

When  s/d  approaches  unity,  we  find  1^(0)  ~  N  1^(0).  One  of  the 

requirements  for  the  multi-element  diode  laser  is  to  develop 

phased  array  designs  in  which  s/d  may  indeed  approach  unity  (see 

the  discussion  in  Section  2.3). 

Having  established  this  simple  relationship,  it  is  possible 

to  use  the  scaling  results  of  Figure  39  to  calculate  how  many 

commercially  available  diodes  are  required  to  increase  the  peak 

2  2 

intensity  from  3.9  kW/cm  in  Reference  21  to  1  MW/cm  ,  thereby 
increasing  the  doubling  efficiency  from  0.15  percent  to  40 
percent,  a  practical  value  typical  of  present  hardware.  The 
output  characteristics  of  the  diode  used  in  Reference  21  are 
shown  in  Table  I ,  along  with  those  of  representative 
commercial  diodes.  Even  allowing  for  transmission  losses  of  as 
much  as  20  percent  in  the  beam  shaping  optics,  it  is  clear  that 
approximately  20  times  more  power  is  available,  within  a  narrower 
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Table  I.  Diode  Laser  Output  Characteristics 


Type 

Unspecified 

RCA 

SG2012 

Laser  Diode 
LD-67 

Wavelength,  nm 

860.0 

904.0 

904.0 

Bandwidth,  nm 

1.2 

3.5 

3.5 

Half-Angle  Beam  Spread  - 

Parallel  to  Junction,  deg 

10.0 

7.5 

7.0 

Normal  to  Junction,  deg 

18.0 

9.0 

10.0 

Peak  Power,  watt 

0. 78a 

20. 0b 

20. 0b 

Pulse  Length,  nsec 

10.0 

200.0 

200.0 

Duty  Cycle,  Percent 

- 

0.1 

0.1 

Reference 

21 

28 

29 

Power  focused  into  nonlinear  crystal 

^Minimum  total  peak  power 

divergence  angle,  than  was  used  in  Reference  21.  Thirteen  of 
these  higher  power  diodes  would  produce  260  W  peak  power  (with 
s/d  =  1).  With  the  same  f-number  as  Reference  21  (i.e.,  a  focal 
length  of  ^1  m)  the  required  peak  intensity  of  1  MW/cm2  could  be 
generated. 

In  summary,  phase  matching  requirements  lead  to  the 

selection  of  KNbO^  as  the  best  nonlinear  crystal.  Using 

established  doubling  performance  data,  a  peak  intensity  of  1  MW/ 
2 

cm  is  projected  to  produce  a  practical  40  percent  doubling 
efficiency.  The  required  peak  intensity  can  be  provided  by  the 
coherent  operation  of  only  thirteen  commercial  diode  lasers. 
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Greater  doubling  efficiency  may  be  achieved  using  intracavity 
doubling.  However,  to  predict  performance  of  that  configuration, 
more  detailed  knowledge  is  required  concerning  such  parameters 
as  single-pass  gain  and  output  coupling. 


4.4  FREQUENCY  DOUBLER  SCALING  ISSUES 

Historically,  the  design  of  high  average  power  frequency 
doubled  systems  has  been  difficult  because  of  the  lack  of  a  full 
understanding  and  appreciation  of  the  physical  mechanisms  and 
phenomena  operating  in  the  nonlinear  crystal  during  high  power 
operation.  This  has  been  aggravated  by  a  paucity  of  reliable 
data  on  the  crystal  properties,  particularly  radiative  absorption, 
that  parameterize  the  power  scaling  relationships.  The  absorp¬ 
tion  measurement  is  difficult  because  the  small  absorptions 
involved  cannot  be  measured  by  conventional  photometric  methods; 
also,  nonlinear  absorption  effects  sometime  arise.  Compounding 
these  difficulties  is  the  variability  in  crystal  properties  with 
vendor . 

The  difficulties  that  have  been  encountered  in  frequency 
doubling  technology  may  be  divided  into  four  broad  areas: 

1.  Efficiency  of  the  doubling  process  and  how  it  is 
limited  by  crystal  type  and  uniformity,  crystal 
environment,  and  laser  beam  properties. 

2.  Measurement  of  the  crystal  properties,  most  importantly 
absorption,  that  parameterize  high  power  operation. 
Absorption  in  the  crystal  is  the  origin  of  the  thermal 
effects  that  limit  output  power. 

3.  Understanding  the  thermal  effects.  This  includes  both 
transient  and  steady-state  effects.  Absorption 
induces  temperature  gradients  in  the  crystal  and  the 
oven  housing  it.  Gradients  outside  the  beam  inter¬ 
action  region  cause  transient  effects  and  instabilit¬ 
ies  in  the  temperature  control.  Gradients  in  the 
interaction  region  degrade  doubling  efficiency  as  they 
disrupt  the  phase  match  conditions.  They  also  cause 
birefringence  in  the  crystal. 
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4. 


Damage  and  deterioration  mechanisms  of  the  materials 
involved.  For  lower  average  power  operation,  surface 
and  bulk  damage  from  high  peak  power  must  be  con¬ 
sidered.  At  high  average  power,  crystal  fracture 
resulting  from  thermal  stress  or  shock  must  also  be 
eliminated.  For  long  life,  material  sensitivity  to 
environmental  factors  must  be  understood  and  controlled. 

Each  of  these  difficulties  has  been  addressed  by  Hughes  in 
earlier  frequency  doubler  designs.  The  important  results  are 
summarized  here  for  each  issue: 

1.  Non-critically  phase  matched  (90  degree  phase  matched) 
materials  have  a  great  advantage  over  angle  phase 
matched  in  that  they  do  not  impose  stringent  beam 
divergence  or  alignment  requirements  on  the  pump  laser 
source  for  high  efficiency.  Higher  doubling  effici¬ 
encies  are  also  obtainable. 

2.  Several  parameters  are  important  for  quantifying  fre¬ 
quency  doubler  power  handling.  Absorption  is  certainly 
the  most  important.  Sensitive  calorimetric  methods  must 
be  utilized  for  measuring  the  small  absorptions  involved, 
so  that  average  power  limits  of  materials  can  be  modelled. 

3.  Thermal  effects  in  the  crystals  are  the  primary  cause  of 
instability  and  limit  the  power  handling.  Transient 
problems  can  be  overcome  by  proper  crystal  oven  design 
and  the  use  of  electro-optic  tuning30  (EOT)  which  pro¬ 
vides  fast  response  to  perturbations  and  always  main¬ 
tains  the  oven  at  the  correct  temperature.  Also,  beam 
shaping31  has  been  shown  to  eliminate  the  average  power 
limit  and  fracture  problems  caused  by  temperature 
gradients  in  the  interaction  region. 

4.  Damage  has  been  found  to  be  a  limiting  factor  in  many 
frequency  doubling  systems.  Many  such  problems  are 
eliminated  by  careful  control  of  contaminants,  particularly 
in  the  crystal  oven  where  outgassing  aggravates  the 
problem.  In  addition,  conservative  design  guidelines  on 
acceptable  flux  levels  have  been  developed  using 
empirical  data. 
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Because  of  their  essential  role  in  the  high  average  power 
frequency  doubler  design,  the  techniques  of  beam  shaping  and  electro¬ 
optic  tuning  will  be  discussed  briefly. 

4.4.1  Beam  Shaping 

The  maximum  output  power  of  any  frequency-doubling  crystal 
is  ultimately  determined  by  the  temperature  gradients  arising 
from  radiative  absorption  in  the  crystal.  This  limit  manifests 
itself  in  two  ways.  First,  because  the  crystal  refractive 
indices  are,  in  general,  temperature  dependent  it  is  impossible 
to  maintain  phase  matching  across  the  entire  beam  profile  in  the 
presence  of  severe  temperature  gradients.  Second,  with  a  con¬ 
ventional  focusing  geometry  which  produces  a  circular  laser  spot 
on  the  nonlinear  crystal,  thermally  induced  stresses  arise  which 
are  essentially  symmetric  about  the  laser  axis.  These  stresses 
depolarize  the  radiation  passing  through  the  nonlinear  crystal; 
the  depolarized  component  is  no  longer  phase  matched  and  is 
effectively  removed  from  the  doubling  process. 

The  solution  to  this  problem  is  called  beam  shaping: 
cylindrical  lenses  shape  the  pump  beam  into  an  ellipse  on  the 
face  of  the  nonlinear  crystal.  A  careful,  but  straightforward 
thermal  analysis31  of  this  configuration  indicates  that  the  maxi¬ 
mum  temperature  difference  across  the  laser  profile  is  reduced 
by  the  factor  (h/w)  relative  to  that  for  a  circular  spot,  where 
h  and  w  are  the  height  and  width  of  the  laser  spot,  respectively. 

For  example,  simply  reducing  h  and  increasing  w  each  by  a  factor 
of  four  would  keep  the  illuminated  area  and  the  intensity  constant, 
but  would  reduce  the  temperature  difference  and  increase  the 
maximum  output  power  by  a  factor  of  sixteen  relative  to  that  for 
a  circular  spot  (h/w  =  1) .  In  addition,  to  the  extent  that  heat 
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flow  in  the  crystal  is  one-dimensional  (using  a  low,  wide  crystal 
geometry  slightly  larger  than  h  and  w  assures  this) ,  a  judicious 
choice  of  pump  beam  polarization  either  a  parallel  or  perpendi¬ 
cular  to  the  direction  of  heat  flow  eliminates  the  depolarization 
problem  when  90°  phase  matching  is  used. 

4.4.2  Electro-Optic  Tuning 

Electro-optic  tuning30  (EOT)  is  a  technique  to  provide  trans¬ 
ient  control  of  frequency  doubling.  It  is  useful  when  long  warmup 
periods  are  undesirable,  or  when  a  high  degree  of  output  amplitude 
stability  is  required.  EOT  is  based  on  the  fact  that  the 
refractive  indices  of  a  nonlinear  crystal,  and  thus  the  phase 
match  condition,  may  be  controlled  much  more  rapidly  by  changing 
the  voltage  across  the  crystal  than  by  changing  its  temperature. 

In  operation,  it  is  initially  assumed  that  the  crystal  is  near, 
but  not  necessarily  exactly  at,  the  optimum  phase  match  tempera¬ 
ture.  A  voltage  of  plus  (minus)  AV  is  applied  to  the  crystal 
for  all  the  odd  (even)  pulses.  Using  a  photodiode  to  measure  the 
second  harmonic  output  energy,  simple  electronic  logic  rapidly 
adjusts  the  average  crystal  voltage  Vq  until  all  output  pulses 
are  of  equal  energy,  an  indication  that  optimum  phase  matching 
has  been  achieved  at  the  peak  of  the  doubling  efficiency  versus 
temperature  or  voltage  curve.  Longer  term  logic  changes  the 

oven  temperature  in  such  a  direction  that  V  never  exceeds  a  few 

o 

kilovolts.  EOT  can  correct  for  temperature  excursions  as  much 
as  a  few  degrees  Kelvin. 
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SECTION  5 


CONCLUSIONS 


This  research  program  was  aimed  at  investigating  several 
technical  areas  related  to  scaling  the  output  power  of  diode 
lasers  by  coherent  coupling  techniques.  We  have  seen  that  the 
far-field  intensity  distribution  of  a  linear  array  of  phase- 
locked  diode  lasers  is  highly  sensitive  to  the  number  spacing 
and  relative  phase  shift  of  adjacent  emitters.  Comparison  with 
recently  published  experimental  results  indicates  that  a  phase 
shift  of  180°  occurs  most  often,  producing  a  slightly  undesir¬ 
able  two-lobed  far-field  pattern.  Computer  simulations  of 
locking  multiple  diode  elements  in  a  phase  conjugate  resonator 
show  that  no  such  problem  arises  in  this  case,  and  that  the 
output  beam  quality  is  determined  solely  by  the  quality  of  the 
resonator  output  mirror.  This  positive  result  is  the  motivation 
for  considering  such  a  coupling  scheme.  Further  computer  simu¬ 
lations  including  practical  features  such  as  waveguiding  effects 
are  recommended  as  a  point  of  further  study.  While  such  calcu¬ 
lations  were  beyond  the  scope  of  the  present  program,  a  first 
step  was  undertaken  in  calculating  the  lowest  order  normal  modes 
of  a  multiple-element  planar  waveguide.  We  have  also  considered 
the  deleterious  effects  on  beam  quality  of  relatively  large  gaps 
between  emitters  (to  accomodate  cooling,  for  example),  and  have 
proposed  several  beam  compacting  schemes  to  minimize  these 
effects.  The  proposed  schemes  are  considerably  more  practical 
in  a  phase  conjugate  resonantor  than  they  would  be  using  other 
coupling  schemes  (e.g.,  evanescent-wave  coupling). 

The  phase  conjugate  mirror  involved  in  constructing  a  phase 
conjugate  resonator  always  involves  some  trade-offs.  While 
devices  based  on  degenerate  four-wave  mixing  (DFWM)  are  concep¬ 
tually  quite  simple,  operate  at  very  low  powers,  and  are  avail¬ 
able  over  a  wide  spectral  range,  their  use  in  the  present  appli¬ 
cation  is  somewhat  undesirable  because  of  a  practical  reason: 
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the  pump  lasers  required  to  operate  the  DFWM  phase  conjugate 
mirror  must  be  diffraction-limited  beams  with  a  total  power 
close  to  that  of  the  multiple-diode  phase  conjugate  resonator. 

In  other  words,  the  pump  laser  must  be  essentially  the  same  as 
the  laser  we  are  trying  to  construct.  This  recognition  leads 
one  to  consider  alternative  phase  conjugation  techniques;  spe¬ 
cifically,  stimulated  Brillouin  scattering  (SBS).  While  SBS 
does  not  require  any  external  pump  laser,  it  does  impart  a  small 
frequency  shift  (<  1  GHz)  on  the  reflected  beam.  Hence,  SBS 
would  be  most  practical  in  pulsed  applications  in  which  the 
accumulated  frequency  shift  could  be  less  than  the  ~  4  nm  laser 
gain  bandwidth  (a  10-cm  long  resonator  could  accomodate  pulses 
of  at  least  1  y sec  duration).  We  have  analyzed  two  representa¬ 
tive  SBS  media  for  potential  application  to  a  diode-based  phase 
conjugate  resonator:  a  fused  silica  fiber,  and  a  CS2~filled 
light  guide.  In  either  case,  we  find  threshold  intensities  of 
~  10-4  -  10“5  w/cm2,  values  which  could  likely  be  achieved  with 
a  high  power  diode  array.  This  fact  supports  the  feasibility  of 
a  diode-based  phase  conjugate  resonator. 

We  have  also  considered  the  practical  issues  involved  in 
achieving  efficient  frequency  doubling  with  a  high  power  diode 
laser.  Non— critically  phase— matched  materials  are  preferred 
because  of  their  wide  angular  tolerance  and  long  interaction 
lengths.  The  baseline  material  selected  is  KNb03 ;  based  on 
published  data  on  the  doubling  efficiency  of  this  material,  we 
find  that  13  coherently  coupled  commercial  diode  lasers  could 
Provide  sufficient  brightness  to  achieve  40%  extracavity 
doubling  efficiency.  Hence,  wavelength  conversion  of  a  high 
power  diode  laser  array  also  seems  practical. 

The  conclusion  is  that  the  basic  concept  of  coupling  multi¬ 
ple  diode  lasers  in  a  phase  conjugate  resonator  appears  feasi¬ 
ble.  Further  theoretical  studies  would  be  appropriate  before 
attempting  an  experimental  demonstation.  These  studies  should 
include  waveguiding  effects,  as  well  as  realistic  performance 


78 


modeling  for  the  phase  conjugate  mirror.  The  analysis  presented 
in  the  body  of  this  report  shows  that  sufficient  brightness  can 
be  obtained  to  allow  efficient  wavelength  conversion  by  nonlin¬ 
ear  processes.  Using  demonstrated  frequency  doubler  power 
scaling  techniques,  an  efficient,  high  power  diode-based  source 
could  produce  radiation  in  the  blue-green  spectral  region  for  a 
variety  of  applications. 


/ 
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